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<N ; Abstract 

The classical Lagrangian of chromodynamics, its quantization in the perturbation 
theory framework, and renormalization form the subject of these lectures. Symme- 
tries of the theory are discussed. The dependence of the coupling constant a s on the 
renormalization scale fi is considered in detail. 

s\ 

1 Introduction 

Oh; 

Many textbooks are exclusively p~H6| or in part [THTT] devoted to quantum chromodynam- 
ics. Quantization of gauge fields is discussed in [12] in detail. Here we'll follow notation 
of [TT]; many calculational details omitted here can be found in this book. References to 
original papers will not be given, except a few cases when materials from such papers was 
directly used in the lectures. 

Quantum chromodynamics (QCD) describes quarks and their interactions. Hadrons 
are bound states of quarks and antiquarks rather than truly elementary particles. Quarks 
have a quantum number called color. We'll present formulas for an arbitrary number of 
C\J ' colors N c ; in the Nature, N c = 3. 
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2 Classical QCD Lagrangian 
2.1 Color group SU{N C ) 

The quark field q l has a color index i G [l,iV c ]. The theory is symmetric with respect to 
transformations 

q l ^>-U t j q j or q ->■ Uq , (1) 
where the matrix U is unitary and has determinant 1: 

U + U = 1 , det U = 1 . (2) 

Such matrices form the group SU(N C ). Quark fields transform according to the funda- 
mental representation of this group. The conjugated quark field q~i = (<? l ) + 7° transforms 
according to the conjugated fundamental representation 

qt -> Wqj or q -> qU + , where W = (IPj)* • (3) 
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The product qq' is invariant with respect to color rotations: 

qq' ->■ qU + Uq' = qq' . (4) 

In other words, Sj is an invariant tensor, its components have the same values (1 or 0) in 
any basis: 

^^rvr/ (■',(■/ 6). (5) 

This tensor describes the color structure of a meson. 

The product of three quark fields SijkqlqiQs = 3) is also invariant: 

tijkqiqlq* -»■ ^ijkU" l v U 3 ' f U h Wi'^'?*' = detU ■ e i/fk >q(q J 2 q% = £tjkQi4^ ■ ( 6 ) 

Here e^k is the unit antisymmetric tensoJE In other words, e^k is an invariant tensor: 

£ijk Si'j'k'Ui 1 Uj-i Uk k = det U + ■ e^k = %fc . (7) 

It describes the color structure of a baryon. The operator with the quantum numbers of 
an antibaryon has the form 

e ijk quq2jhk -»> e ijk q u q 2j q 3k ■ (8) 

I.e., e^ k is also an invariant tensor. 

The matrix of an infinitesimal color rotation has the form 

U = 1 + ia a t a , (9) 

where a a are infinitesimal parameters, and the matrices t a are called the generators of the 
fundamental representation of the group SU(N C ). The properties ([2]) of matrices U imply 
that the generators are hermitian and traceless: 

U+U = 1 +ia a (t a - (t a ) + ) = 1 (t a ) + = t\ , , 

det 17 = 1 + ia a Trt a = 1 TrT = 0. ^ 

The trace 

Tr t a t b = T F 6 a \ (11) 

where Tp is a normalization constant (usually Tp = | is chosen, but we'll write formulas 
with an arbitrary Tp). 

How many linearly independent traceless hermitian matrices t a exist? The space of 
hermitian iV c x N c matrices has dimension N%; vanishing of the trace is one additional 
condition. Therefore, the number of the generators t a , which form a basis in the space of 
traceless hermitian matrices, is equal to N% — 1. 



1 In the case of N c colors it has N c indices, and the invariant product contains N c quark fields; a baryon 
consists of N c quarks. 
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The commutator [t a ,£ fe ] is antihermitian and traceless, and hence 

[t\t h ] = if abc t c , (12) 

where 

/^ = _L Tr[0 y (13) 

are called the structure constants of the group SU(N C ). 

Let's consider the quantities A a = qt a q' . They transform under color rotations as 

A a -> qU + t a Uq' = U ab A b , (14) 

where 

U + t a U = U ab t b , (15) 

and hence 

V ab = 1 TrU + t a m b_ ^ 

Tf 

The quantities A a (there are N% — 1 of them) transform according to a representation of 
the group SU(N C ); it is called the adjoint representation. 

Components of the generators {t a ) l j have identical values in any basis: 

(ifYj -> U ab U i ll U/{t b ) i ' r = {t^j , (17) 

hence they can be regarded an invariant tensor. 

The quantities A a transform under infinitesimal color rotations as 

where 

U ab = S ab + ia c {t c ) ab , (19) 
and the generators of the adjoint representation are 

{t c ) ab = if acb . (20) 

Generators of any representation must satisfy the commutation relation ffl2l . In par- 
ticular, the relation 

(t a ) dc (t b ) ce - (t b ) dc (t a ) ce = if abc (t c ) de . (21) 

must hold for the generators (1201) of the adjoint representation. It can be easily derived 
from the Jacobi identity 

[t a , [t b , t d }} + [t b , [t d , t a }} + [t d , [t a , t b ]} = (22) 

(expand all commutators, and all terms will cancel). Expressing all commutators in the 
left-hand side of according to the formula ffTi2"|) . we obtain 

(if bdc if ace + if^if ce + if abc if dce ) t e = 0, (23) 

and hence (12T|) follows. 
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2.2 Local color symmetry and the QCD Lagrangian 

The free quark field Lagrangian 

L = q{iYdf, - m)q (24) 

is invariant with respect to global color rotations q(x) — > Uq(x) (where the matrix U does 
not depend on x). How to make it invariant with respect to local (gauge) transformations 
q(x) — > U(x)q(x)l To this end, the ordinary derivative d^q should be replaced by the 
covariant one D^q: 

D fl q = (d»-igA fl )q, A, = Ap a . (25) 

Here A a ^(x) is the gluon field, and g is the coupling constant. When the quark field 
transforms as q — > q' = Uq, the gluon one transforms too: A^ — > A'^. This transformation 
should be constructed in such a way that D^q transforms in the same way as q: D^q — >• 
D'^q' = UD^q. Therefore, 

(d li -igA' li )Uq = U(d li -igA li )q, 
or dfJJ — igA'^U = —igUA^. We arrive at the transformation law of the gluon field 

4 = UA^W 1 - '-(d^U- 1 . (26) 

Infinitesimal transformations of the quark and gluon fields have the form 

q(x) q\x) = (1 + ia a (x)t a )q(x) , 
A;(x) A»(x) = Al(x) + - g Dfa\x) , (2?) 

where the covariant derivative acting on an object in the adjoint representation is 

Df = 5 ab d^ - ig(t c ) ab A c ^ . (28) 
The expression D v \q also transforms as q: [D' , D' v \q' = U[D^, D v \q. Let's calculate 

it: 

[Dp, D v \q = d^q - ig(d lx A v )q - igA^d^q - igA^d u q - g 2 A ll A v q 
- d v d^q + ig(d u A fJ )q + igA^d v q + igA v d^q + g 2 A v A tl q . 

All derivatives have canceled, and the result is —igG^ u q where 

Gfxv = d^A v — dyA^ — ig[A^, A v \ = G^ u t , 
G% = d,Al - d u Al + 9 r hc A\Al 

is called the gluon field strength. It transforms in a simple way 

UG, u U~ l , G% U ab G% ; (30) 



(29) 



there is no additive term here, in contrast to (1261) . 

Now, at last, we are ready to write down the complete QCD Lagrangian. It contains 
rif kinds (flavors) of quark fields qj and the gluon field: 

L = L q + L A . (31) 

The first term describes free quark fields and their interaction with gluons: 

/ 

The second term is the gluon field Lagrangian: 

L A = -^-Tr G^CT = ~G%G^ . (33) 

It is gauge invariant due to fl30l) . In contrast to the photon field Lagrangian in QED, it 
contains, in addition to terms quadratic in A^, also cubic and quartic terms. The gluon 
field is non-linear, it interacts with itself. 



2.3 Symmetries 

The QCD Lagrangian is symmetric with respect to translations and Lorentz transforma- 
tions, as well as discrete transformations^ P, C, T. It is also symmetric with respect to 
local (gauge) color transformations. 

The QCD Lagrangian is symmetric with respect to phase rotations of all quark fields: 

q f ->• e ia q f « (1 + ia)q f . (34) 

This U(l) symmetry leads to conservation of the total number of quarks minus antiquarks, 
i. e. of the baryon charge. 

If several kinds (flavors) of quarks have equal masses (m/ = m), a wider symmetry 
appears: 

Qf -> Uff'Vf > (35) 

where U is an arbitrary unitary matrix [U + U = 1). Any such matrix can be written as 
U = e ia Uo where det Uq = 1. In other words, the group of unitary transformations is a 
direct product U{nf) = U(l) x SU(nf). Infinitesimal transformations have the form 

U = 1 + ia + ia a r a , (36) 

2 So called tf-term G^G atlv (where G« u = \e^ al3 G aaP ) is not symmetric with respect to P (and CP). 
However, it is the full divergence of a (non gauge-invariant) axial vector. Therefore adding this term (with 
some coefficient) to the Lagrangian changes nothing in classical theory. In quantum theory, the $-term is 
inessential in perturbation theory, but it changes the behavior of the theory due to nonperturbative effects, 
leading to P (and CP) violation in QCD. Such violations have not been seen experimentally; therefore we 
shall not discuss the i?-term. 
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where the SU(jif) generators are hermitian matrices satisfying the condition Trr a = 0. 

In the Nature, masses of various quark flavors are not particularly close to each other. 
However, u and d quarks have masses much smaller than the characteristic QCD energy 
scale (we'll discuss this scale later). Both of these masses can be neglected to a good 
accuracy, and the 577(2) symmetry (called isospin) has a good accuracy (of order 1%). 
The s quark mass is smaller than the characteristic QCD scale but not much so, and the 
flavor SU(3) symmetry has substantially lower accuracy. 

Under a stronger assumption that masses of several quark flavors rrif — 0, left and right 
quarks 

1 ± 75 

Qf = QLf + QRf , Ql,r = — g — q ' 1s<1l,r = ±Ql,r (37) 
live their own lives without transforming to each other: 

L q = QLfi'fD ll q Lf + QRfil^D^Rf (38) 
/ / 

(the mass term mqq = m(qi J q R J r qRqi) transforms left quarks to right ones and vice versa). 
The theory has a larger symmetry U{nf)i x U(nf) R ; its infinitesimal transformations are 

q L -> (1 + ia L + ia a L r a )q L , q R ; -> (1 + ia R + ia a R r a )q L . (39) 

They can be re-written as 

q -> (1 + ia v + iayT a + ia A ^ 5 + ia A r a ^f 5 )q , (40) 

this corresponds to the U(rif)y x U{nf)A symmetry. As already discussed, this symmetry 
is quite good for u and d quarks (n/ = 2), and substantially less accurate if s quark is 
added (n/ = 3). 

If all quarks are massless, then the Lagrangian contains no dimensional parameters, 
and it is symmetric with respect to scale transformations 

x^Xx^, A^\- l A^ q^\- 3/2 q. (41) 

For a wide class of field theories one can prove that the scale invariance implies invariance 
with respect to inversion 

x» -> % (42) 

x A 

(in an infinitesimal neighborhood of each point it is a scale transformation). Performing 
inversion, then translation by a, and then again inversion produces a special conformal 
transformation 

x M + a^x 2 , , 

x M -)■ — . (43) 

1 + 2a ■ x + a 2 x 2 

These transformations, together with scale ones, translations, and Lorentz transformations 
form the conformal group. The classical massless QCD is invariant with respect to this 
group. 

Not all symmetries of the classical theory survive in the quantum one (Tabled]). 
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Table 1: Symmetries of massless QCD in classical and quantum theory. 



Group 


Classical theory 


Quantum theory 


translations 
Lorentz 




conformal 




anomaly 


SU(N C ) 


local 


U(l) 
SU(n f ) 

U(1)a 
SU(n f ) A 






anomaly 

spontaneously broken 


P 
C 
T 


discrete 



3 Quantization 

3.1 Faddeev— Popov ghosts 

It is convenient to use the functional integration method to quantize gauge theories. The 
correlator of two operators 0(x) and 0(y) (we assume them to be gauge invariant) is 
written as 

fU,n,,dA c " l (x)e i -f Ld ' lx O(x)0(y) 1 1 d 2 Z\i] 
<T{0{x) ) 0{y)}>= J - llx ' a ^ ^ ' — " 



, (44) 

3=0 



IU x , a ,,dA^x) e*f Ld4 * ~ Z[j] 5j(x) 5j(y) 
where the generation functional 

m = J U dA >^ elIiL+1 ° )d4x - ( 45 ) 



x,a,fi 



To make formulas shorter, we'll consider gluodynamics (QCD without quarks); including 
quark fields introduces no extra difficulties, one just have to add integration in fermionic 
(anticommuting) fields. 

In the case of gauge fields a problem appears: a single physical field configuration is 
taken into account infinitely many times in the integral. All potentials obtained from a 
given one by gauge transformations A — > A u form an orbit of the gauge group; physically, 
they describe a single field configuration. It would be nice to include it in the functional 
integral just oncqj To this end, one has to fix a gauge - - to require some conditions 
Q a (A u (x)) = For any A(x) this equation should have a unique solution U(x). I.e., the 
"surface" G = should intersect each orbit of the gauge group at a single point (Fig. [T|). 



3 This is necessary in perturbation theory to define the gluon propagator. Gauge fixing is not needed 
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Figure 1: Gauge fixing. 



For example, the Lorenz gaugcl 

G a {A{x)) = VAfa) . (46) 

is often used. The axial gauge G a (A(x)) = n^A a ^(x) (with some fixed vector n) and the 
fixed-point (Fock-Schwinger) gauge G a (A(x)) = x^A®(x) are also popular. 
Let's define the Faddeev-Popov determinant A [A] by the formula 

A" 1 [4 = I \{dU{x) l[6(G a (A u (x))), (47) 

x x,a 

where dU is the invariant integration measure on the group (it satisfies the condition 
d(UoU) = dU); for infinitesimal transformations dU = fT da a . Near the surface G(A U °) = 
0, variations of G at infinitesimal gauge transformations are linear in their parameters: 



so that 



8G{A(x)) = Ma{x) , (48) 
A~ l [A] = J Y[da{x)5{Ma{x)) = 1/detM, (49) 



i. e, A[A] is the determinant of the operator M. For example, for the Lorenz gauge 
G a (A(x)) = we obtain from P7jl 

SG a (x) = -dPD*a\x) => M = ~d»Df] (50) 
9 9 



in some approaches, e.g., the lattice QCD. QCD in Euclidean space-time (obtained by the substitution 
t = —He) is considered. The oscillating expiS becomes exp(— Se) where the Euclidean action Se is 
positive. Continuous space-time is replaced by a discrete 4-dimensional lattice, the exact gauge invariance 
is preserved. Random field configurations are generated with probability exp(— Se)', fields belonging to an 
orbit of the gauge group are generated equiprobably. 

4 The solution U of the equation G(A U ) = for this gauge is not unique (Gribov copies), if one considers 
gauge transformations sufficiently far from the identical one. This problem is not essential for construction 
of perturbation theory, because it is sufficient to consider infinitesimal gauge transformations. 
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for the axial gauge 

1 () ab 

M = -n»D* = —n"d, , (51) 
9 * 9 



invariant: 



because n^A^ = due to the gauge condition. The Faddeev-Popov determinant is gauge 
A" 1 ^ ] = [Y[dU l[8(G a (A UoU (x))) 

x x,a 

]JD(U U) l[6(Cr(A u ° u (x))) = A-'IA] . 

x x,a 

Now we insert the unit factor (1471) in the integrand (145D : 
Z[J}= I Y[dA(x) e lS[A] 

J X 

]JdU(x) ]JdA(x)A[A} ]j5(G(A u (x)))e iS[A] 
Yl I dU J x f \\dA{x)A[A] ]j5(G(A(x))) e lS[A] . (52) 

X / X X 

Here the first factor is an (infinite) constant, it cancels in the ratio ( |44"|) and may be omitted. 
We have arrived at the functional integral in a fixed gauge. 

The only integral which any physicist can calculate (even if awakened in the middle of 
night) is the Gaussian one: 

f dz*dze- az * z ~ -. 
J a 

The result is obvious by dimensionality. In the multidimensional case the determinant 
appears because the matrix M can be diagonalized: 

U^dzte-W-^. (53) 

i 

Integration in a fermion (anticommuting) variable c is defined as 

dc = 0, cdc=l (54) 



(hence fermion variables are always dimensionless) . From c 2 = 

e -ac*c = ! _ flc * C) 

so that 

dc* dce~ ac * c = a . 



In the multidimensional case 



I Y[ dc* da e~ M ^ ~ det M . (55) 

The functional integral ( 1521) is inconvenient because it contains A[A] = detM. It can 
be easily written as an integral in an auxiliary fermion field: 

A [A] = det M = f]J de ( x ) d(f ( x ) e ifLcd4x , L c = -c a M ab c b . (56) 

x,a 

The scalar fermion field c a (belonging to the adjoint representation of the color group, 
just like the gluon) is called the Faddeev-Popov ghost field, and c a is the antighost field. 
Antighosts are conventionally considered to be antiparticles of ghosts, though c a and c a 
often appear in formulas in non-symmetric ways. 

In the axial gauge (15T|) A [A] = det M does not depend on A, and this constant factor 
may be omitted; there is no need to introduce ghosts. They can be introduced, of course, 
but the Lagrangian (J56l) shows that they don't interact with gluons, and thus influence 
nothing. The same is true for the fixed-point gauge. 

In the generalized Lorenz gauge G a (A(x)) = d^A a ^(x) — u a (x) we have, up to an inessen- 
tial constant factor, M = d^D^, and therefore 

L c = -c a d»Dfc b => {d»c a )Dfc b . (57) 

A full derivative has been omitted in the last form; note that the ghost and antighost fields 
appear non-symmetrically — the derivative of c is covariant while that of c is the ordinary 
one. 

The generating functional 

Z[J] = / IJ dAa ^ d ° a ^ II 5 ^ A K x ) - ua ( x )) e i ^ LA+Lc+JO)dix (58) 

x,a x,a 

is gauge invariant; in particular, it does not depend on u a (x). Let's integrate it in 
Tlx a doJ a (x) with the weight exp [— ^ / co a (x)co a (x) d 4 x\ : 

Z[J\ = I Y[ dA a (x) dc a (x) dc a (x) e l ^ L+J °^ x , (59) 

x,a 

where the QCD Lagrangian (without quarks) in the covariant gauge L = L A + L a + L c 
contains 3 terms: the gluon field Lagrangian La', the gauge- fixing term L a ; and the ghost 
field Lagrangian L c , 

L A = -\G%G a ^ , L a = -1 (d^Atf , L c = {d»c a )D ab c b . (60) 
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If there are quarks, their Lagrangian L q ( )32|) should be added, as well as extra integrations 
in quark fields. In quantum electrodynamics ghosts don't interact with photons, and hence 
can be ignored. 

The Lagrangian ( 1601) obtained as a result of gauge fixing is, naturally, not gauge invari- 
ant. However, a trace of gauge invariance is left: it is invariant with respect to transfor- 
mations 

6A% = X + D ab c b , 5c a = —X + d^A a u , 5c a = -?-f abc X + c b c c , (61) 



where A is an anticommuting (fermion) parameter. This supersymmetry (relating boson 
and fermion fields) is called the BRST symmetry. 

3.2 Feynman rules 

The quark propagator has the usual form 

. >- • =iS (p), S (p) = = 4^- 2 , (62) 

p p — m p z — m z 

where the unit color matrix (in the fundamental representation) is assumed. 

It is not possible to obtain the gluon propagator from the quadratic part of the La- 
grangian La'- the matrix which should be inverted is not invertible. The gauge fixing 
procedure is needed to overcome this problem. In the covariant gauge (1601) the quadratic 
part of La + L a gives the gluon propagator 



a 



,,«v/wwww^ -iS ab D° (v) D° (v) = — 

\l p v iiv \r) i ^ iivKV) ^2 



g^u - (i - do) — 5- 

p2 



(63) 



The ghost propagator 



- -> - - - 6 = it ab G {p) , G (p) = - 2 , (64) 
p p z 



as well as the gluon one (1631) . has the color structure 5 ab — the unit matrix in the adjoint 
representation. 

The quark-gluon vertex (see ( 152]) ) 

a 



> < > — = t a x ig ^ (65) 

has the color structure t a ; otherwise it has the same form as the electron-photon vertex in 
quantum electrodynamics. 
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The gluon field Lagrangian La (|33|) contains, in addition to quadratic terms, also ones 
cubic and quartic in A. They produce three- and four-gluon vertices. The three-gluon 
vertex has the form 





= i/™ xig V^ 3 (p 1 ,p 2 ,p 3 ), 



(66) 



It is written as the product of the color structure if' 



and the tensor structure. To 



do such a factorization, one has to choose a "rotation direction" around the three-gluon 
vertex (clockwise in the formula (166]) ) which determines the order of the color indices in 
j-aia 2 a 3 ag we vj ag order of the indices and the momenta in V tlltJ-2f13 (pi, p 2 , p%) ■ Inverting 
this "rotation direction" changes the signs of both the color structure and the tensor one. 
The three-gluon vertex does remains unchanged — it does not depend on an arbitrary 
choice of the "rotation direction" . This choice is required only for factorizing into the color 
structure and the tensor one; it is essential that their "rotation directions" coincide. 

The four-gluon vertex does not factorize into the color structure and the tensor one - 
it contains terms with three different color structures. This does not allow one to separate 
calculation of a diagram into two independent sub-problems — calculation of the color 
factor and of the remaining part of the diagram. This is inconvenient for writing programs 
to automatize such calculations. Therefore, authors of several such programs invented the 
following trick. Let us declare that there is no four-gluon vertex in QCD; instead, there is 
a new particle interacting with gluons: 




The propagator of this particle doesn't depend on p: 



In coordinate space it is proportional to 5(x), i.e., this particle does not propagate, and 
all four gluons interact in one point. Interaction of this particle with gluons has the form 

v 



a%^£b = -5 ab (g^ a g u(3 



V a 2 



~9' 





(69) 
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The sum (1671) correctly reproduces the four-gluon vertex following from the Lagrangian 
La dSSjIE The number of diagrams increases, but each of them is the product of a color 
factor and a "colorless" part. 

Finally, the ghost-gluon vertex has the form 



It contains the momentum of the outgoing ghost but not of the incoming one because of the 
asymmetric form of the Lagrangian fl60|) . In the color structure, the "rotation direction" 
is the incoming ghost — > the outgoing ghost — > the gluon. 

The color factor of any diagram can be calculated using the Cvitanovic algorithm. It 
is described in the textbook [TT] . 

4 Renormalization 



4.1 MS scheme 

Many perturbation-theory diagrams containing loops diverge at large loop momenta (ul- 
traviolet divergences). Because of this, expressions for physical quantities via parameters 
of the Lagrangian make no sense (contain infinite integrals). However, this does not mean 
that the theory is senseless. The requirement is different: expressions for physical quanti- 
ties via other physical quantities must not contain divergences. Re-expressing results of the 
theory (which contain bare parameters of the Lagrangian) via physical (i.e., measurable, 
at least in principle) quantities is called renormalization, and it is physically necessary. 
Intermediate results of perturbation theory, however, contain divergences. In order to give 
them a meaning, it is necessary to introduce a regularization, i. e. to modify the theory 
in such a way that divergences disappear. After re-expressing the result for a physical 
quantity via renormalized parameters one can remove the regularization. 

The choice of regularization is not unique. A good regularization should preserve as 
many symmetries of the theory as possible, because each broken symmetry leads to con- 
siderable complications of intermediate calculations. In many cases (including QCD) it 
happens to be impossible to preserve all symmetries of the classical theory. When a regu- 
larization breaks some symmetry, intermediate calculations are non-symmetric (and hence 
more complicated); after renormalization and removing the regularization, the symmetry 
of the final result is usually restored. However, there exist exceptions. Some symmetries 
are not restored after removing the regularization, they are called anomalous. I.e., these 
symmetries of the classical theory are not symmetries of the quantum theory. 

One can also prove this equivalence using functional integration (see |13|). We remove the terms 
quartic in A from La and introduce an antisymmetric tensor field with the Lagrangian — \t a>iv t a ilv + 
■^f abc t a ^ v A^AI (producing the Feynman rules (1691) ). It is easy to calculate the functional integral 
in this field, and the QCD generating functional with the full La is reproduced. 



[I b 




(70) 



P 
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In the case of gauge theories, including QCD, it is most important to preserve the 
gauge invariance. For example, the lattice regularization used for numerical Monte-Carlo 
calculations preserves it. However, this regularization breaks translational and Lorentz 
invariance (Lorentz symmetry restoration in numerical results is one of the ways to estimate 
systematic errors). Because of this, the lattice regularization is inconvenient for analytical 
calculations in perturbation theory. 

In practice, the most widely used regularization is the dimensional one. The space-time 
dimensionality is considered an arbitrary quantity d = 4 — 2e instead of 4. Removing the 
regularization at the end of calculations means taking the limit e — > 0; intermediate expres- 
sions contain l/e n divergences. Dimensional regularization preserves most symmetries of 
the classical QCD Lagrangian, including the gauge and Lorentz invariance ((^-dimensional). 
However, it breaks the axial symmetries and the scale (and hence conformal) one, which 
are present in the classical QCD Lagrangian with massless quarks. The scale symmetry 
and the flavor-singlet U(1)a symmetry appear to be anomalous, i. e. they are absent in the 
quantum theory. 

Now we'll discuss renormalization of QCD in detail. For simplicity, let all n/ quark 
flavors be massless (quark masses are discussed in Sect. [6]). The Lagrangian is expressed 
via bare fields and bare parameters; in the covariant gauge 

L = J2 QodYD,q 0i - -G^ u G a r - ^ (<W) 2 + (d^)(D,c a ) , (71) 

i 

where 

D^qo = (<9 M - ig Aon) q , A 0fl = A^t" , 
[Dp, D u ]q = —igoG 0fJjU q , G 0fiu = G^^t " , 
Go^v = d^A^ — d u A^ + g f a °A ^Aq 1 , , 
DA = (d,5 ab - igoA$)c b , A* = A c 0/1 {t c ) ab . 

The renormalized fields and parameters are related to the bare ones by renormalization 
constants: 

q = Z 1 q / 2 q, A = Z l J 2 A, a = Z A a, g = Z 1 J 2 g (72) 

(we shall soon see why renormalization of the gluon field and the gauge parameter a is 
determined by a single constant Za)- In the MS scheme renormalization constants have 
the form 

z ^ s) = l+ l^ + ^ + -)^) + --- (73) 

In dimensional regularization the coupling constant g is dimensional (this breaks the scale 
invariance). Indeed, the Lagrangian dimensionality is [L] = d, because the action must be 
dimensionless; hence the fields and go have the dimensionalities [Aq] = 1 — e, [qo] = 3/2 — e, 
[go] = e. In the formula ( 173|) a s must be exactly dimensionless. Therefore we are forced to 
introduce a renormalization scale \i with dimensionality of energy: 

4tt ^ (4vr) d / 2 ' 1 ' 
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The name MS means minimal subtraction: minimal renormalization constants (175)1 contain 
only negative powers of e necessary for removing divergences and don't contain zero and 
positive power^l In practice, the expression for g$ via a s (fi), 

91 -^^-Z a (a s (^, (75) 



(47r) d / 2 ^ 4vr 

is used more often. First we calculate something from Feynman diagrams, results contain 
powers of go] then we re-express results via the renormalized quantity a s (fi). 



4.2 The gluon field 

The gluon propagator has the structure 



+ 



+ 



+ ■ ■ 



+ 



iD%(p) + {-i)Dl a { P )iU^{p){-i)D%{j>) 

)Dl a (p)M a P(p)(-i)D° (p)iW 5 (p)(-i)Dl (p) + 



(76) 



where the gluon self energy i5 ab U^{p) is the sum of all one particle irreducible diagrams 
(which cannot be cut into two disconnected pieces by cutting a single gluon line). This 
series can be re-written as an equation: 



DM = DIM + Dl a {p)^ip)DM ■ 
For each tensor of the form 



(77) 



A 



A 



p 



I p2 



it is convenient to introduce the inverse tensor 



.4 



A- 1 



PyPv 
p2 



+ An 1 



p 2 



satisfying 



Then the equation fITTj) can be re-written in the form 

d; 1 (p) = (d°)-J( p )-ti, u (p). 

In a moment we shall derive the identity 

iV(pK = o 



(7* 



(79) 



6 The bar means the modification of the original MS scheme introducing the exponent with the Euler 
constant 7 and the power d/2 instead of 2 in the denominator — these changes make perturbative formulas 
considerably simpler. 
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which leads to 

n M ,(p) = (pV-^) n (p 2 )- (so) 

Therefore, the gluon propagator has the form 

1 



DM 



p 2 (i - rv)) [ M P 



p^p 

9iw 



,2 



PpVv 

+ a 



(P 



2\2 



There are no corrections to the longitudinal part of the propagator. The renormalized 
propagator (related to the bare one by D^ u {p) = ZA(a{(i))D r ^ u (p; /x)) is equal to 



_ PjJK 



+ a 0*)jjff- ( 82 ) 



The minimal ( 1731) renormalization constant Za{ol) is tuned to make the transverse part of 
the renormalized propagator 

irjj?- } v) = zi\a{ti) 1 



p 2 (l - Tl(p 2 )) 



finite at e — >■ 0. But the longitudinal part of (152"]) (containing a(/i) = Z^ 1 (a(/x))a ) also 
must be finite. This is the reason why renormalization of is determined by the same 
constant Za as that of the gluon field (!72|) . 

In quantum electrodynamics, the property ( 175]) follows from the Ward identities, and 
its proof is very simple (see, e.g., [IT]). In quantum chromodynamics, instead of simple 
Ward identities, more complicated Slavnov-Taylor identities appear; transversality of the 
gluon self energy (I79p follows from the simplest of these identities. Let's start from the 
obvious equality 



<T{d»A;(x),c b (y)}> = 

(single ghosts cannot be produced or disappear, as follows from the Lagrangian (ITT]) ). 
Variation of this equality under the BRST transformation ( 16T]) is 

<T{d^(x), d v A h v {y)}> - a<T{d»D; c c c (x), c\y)}> = . 

Using the equation of motion for the ghost field d^D^c? = 0, we arrive at the Slavnov- 
Taylor identity 

<T{d»A«(x),d»A b v (y)}> = 0. (83) 

The derivative 

±.±. <T{A>) , Alm> 

does not vanish: terms from differentiating the ^-function in the T-product remain. These 
terms contain an equal-time commutator of A^ (x) and A^y); it is fixed by the canonical 
quantization of the gluon field, and thus is the same in the interacting theory and in the 
free one with g = 0: 

P >1 p u D, u (p)=P fl P U D%(p). 
16 



Hence (179|) follows. 

The gluon self energy in the one- loop approximation is given by three diagrams (Fig. [2]). 
The quark loop contribution has the structure ( l80l) . and can be easily obtained from the 
QED result. The gluon and ghost loop contributions taken separately are not transverse; 
however, their sum has the correct structure (ISUl) . Details of the calculation can be found 
in ITT]. The result is 



n(p : 



2, 9 2 o(-pT £ 



(47r) d / 2 2(d-l) 
+ C A 



-AT F n f (d-2) 



3d-2 + (d- l)(2d - 7)£ - -(d - l)(d - 4)£ 2 



where £ = 1 — a , 



Gi 



2gi 



(d — 3)(d — 4) 



5-1 



r(i + g )r 2 (i- £ ) 

r(l-2e) 



*4) 



(85) 



o o o 



Figure 2: The gluon self energy at one loop. 

The transverse part of the gluon propagator (1ST]) , expressed via the renormalized quan- 
tities a s (/i) (1751) and a(/i), is 



P 2 £±(P 2 ) 



e" i£ e 7£ - 



<7i 



16(1 — e)T F nf 



Aire ' 4(l-2e)(3-2e) 
- (e(3 - 2e)a 2 (fi) - 2(3 - 2e)(l - 3e)a(/i) + 26 - 37e + 7e 2 ) CU 

where L = log(— p 2 / fj 2 ). Expanding in e (e ie gi = 1 + 0(e 2 )), we get 



p 2 D ± (p 2 



47r£ 



+ l' 9 " 2 + 1 3 8 6 a + 9 V -^ nj ) e 



The result must have the form p 2 D±(p 2 ) = Za(ol s (h), a(fi))p 2 D r L (p 2 ] //) where D r ± (p 2 ; /i) is 
finite at e — >• 0. Therefore, at one loop 



Z A (a s ,a) = 1 



a s 
Aire 



(86) 
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4.3 Quark fields 

The quark propagator has the structure 



+ 



+ 



+ 



iS(p) = iS {p) + iS (p)(-i)Z(p)iS {p) (87) 
+ iS (p){-i)Z{p)iS {p){-i)Z{p)iS (p) + ■■■ 

where the quark self energy —iL{p) is the sum of all one particle irreducible diagrams 
(which cannot be cut into two disconnected pieces by cutting a single quark line). This 
series can be re-written as an equation 

Sip) = S (p) + S (pMp)S(p) ; (88) 

its solution is 

S(P) = q _ 1( * =r ( : • (89) 
S (P)-S(p) 

For a massless quark, S(p) = p"Ey(p 2 ) from helicity conservation, and 

The quark self energy in the one-loop approximation is given by the diagram in Fig. [3J 
Details of the calculation can be found in [IT], the result is 



CLqGi . 




(91) 



Figure 3: The quark self energy at one loop. 
The quark propagator fl9"Ul) . expressed via the renormalized quantities and expanded in 



e, is 



4?r * (d-3)(d-A) 

= l-C F ^-aQi)e- u (l + e + ■••)• 

It must have the form Z q (a(fi), a(/j,))p'S r (p; p) where S r (p; p) is finite at e — > 0. Therefore, 
at one loop 

ZJa,a) = l-C F a—?-. (92) 
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4.4 The ghost field 

The ghost propagator is 



G(p) 



p 2 — X(j3 2 ) 

The ghost self energy in the one-loop approximation (Fig. [3]) is (see [TT] ) 



s(p 2 ) = - 7 g 



1^ 0o(~P 



It „,2\l-e 



-Gi [<2 — 1 — (d — 3)a ] . 



(93) 



(94) 



4 A (47r) d / 2 

Re-expressing the propagator via the renormalized quantities and expanding in e, we get 

a s (/x) _ Le 3-a + 4e 



and hence 



p 2 G(p) = 1 + C A ^e 



ZJa s , a) = 1 + C 



4:718 4 

3 — a a. 



4 47T£ 




Figure 4: The ghost self energy at one loop. 



(95) 



5 Asymptotic freedom 

In order to obtain the renormalization constant Z a , it is necessary to consider a vertex 
function and propagators of all the fields entering this vertex. It does not matter which 
vertex to choose, because all QCD vertices are determined by a single coupling constant g. 
We shall consider the quark-gluon vertex. This is the sum of all one particle irreducible di- 
agrams (which cannot be separated into two parts by cutting a single line), the propagators 
of the external particles are not included: 




tg t a T»(p,p') , T»(p,p') =Y + A M (P,P') 



(96) 



p p 
(A M starts from one loop). 
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The vertex function expressed via the renormalized quantities should be equal to T M = 
Z-pT^, where Z-p is a minimal ([75]) renormalization constant, and the renormalized vertex 
is finite at e — > 0. 

In order to obtain a scattering amplitude (an element of the S- matrix), one should 
calculate the corresponding vertex function and multiply it by the field renormalization 
constants Z^ 2 for each external particle i. This is called the LSZ reduction formula. We 
shall not derive it; it can be intuitively understood in the following way. In fact, there are 
no external lines, only propagators. Suppose we study photon scattering in the laboratory. 
Even if this photon was emitted in a far star (Fig. [5]), there is a photon propagator from 
the star to the laboratory. The bare photon propagator contains the factor Za- We split 
it into Z 1 / 2 ■ Z 1 / 2 : , and put one factor Z 1 / 2 into the emission process in the far star, and 

1/2 

the other factor Z A into the scattering process in the laboratory. 




Figure 5: Scattering of a photon emitted in a far star. 

The physical matrix element goTZ q Z][ 2 = gT r Zl/ 2 Z-pZ q Z^ 2 must be finite at e — » 0. 
Therefore the product Z]j 2 Z^ZqZ 1 ^ 2 must be finite. But the only minimal ( 1731) renormal- 
ization constant finite at e — > is 1, and hence 

Z a = (ZrZ^Z^ . (97) 

In QED ZrZ q = 1 because of the Ward identities, and it is sufficient to know Za- In QCD 
all three factors are necessary. 

The quark-gluon vertex in the one-loop approximation is given by two diagrams (Fig.E]). 
In order to obtain Zp, it is sufficient to know ultraviolet divergences (l/e parts) of these 
diagrams; they don't depend on external momenta. Details of the calculation can be found 
in [11], the results for these two diagrams are 

and hence 

Z r = l+(c F0 + C,^)^. (98) 
From flSSD and ([92]) we obtain 

z r z q = i + c A ^^- £ . 
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Figure 6: The quark-gluon vertex at one loop. 

The color structure Cp has canceled, in accordance with QED expectations. Finally, 
taking (|86p into account, the renormalization constant Z a f[9"T|) is 

— C A - -T F n f ) . (99) 
3 A 3 V 4tt£ v ' 

It does not depend on the gauge parameter a, this is an important check of the calculation. 
It can be obtained from some other vertex, e.g., the ghost-gluon one (this derivation is 
slightly shorter, see [TTj). 

Dependence of a s (//) on the renormalization scale /i is determined by the renormal- 
ization group equation. The bare coupling constant g\ does not depend on \i. Therefore, 
differentiating the definition (|74p in d log /i, we obtain 

d ^M = -2e-2P(aM), (100) 
d log fj, 

where the /3-function is defined as 

flo.00) = \ d]o ^° [aM) . (ioi) 

2 a log /i 
For a minimal renormalization constant 

Z a (a s ) = 1 + z x ^- H 

we obtain from (11 01 1) with one-loop accuracy 

/?(a s ) = Ay 1 H = -^ly 1 H 

4tt 47T 

This means that the renormalization constant Z a has the form 

Z a (a s ) = 1 - /3 -r^ h — 

47T£ 

From fl99l) we conclude that 

Pv = ^C A -^T F n f . (102) 

For n/ < = y this means that f3(a s ) > at small a s , where perturbation theory 
is applicable. In the Nature rif = 6 (or less if we work at low energies where the existence 
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of heavy quarks can be neglected), so that this regime is realized: a s decreases when the 
characteristic energy scale fi increases (or characteristic distances decrease). This behavior 
is called asymptotic freedom; it is opposite to screening which is observed in QED. There 
the charge decreases when distances increase, i. e. \i decreases. 
The renormalization group equation with e = 0, 

dloga a (fi) 

— = -2p[a a (ji)) , 

d log /i 

can be easily solved if the one-loop approximation for /3(a s ) is used: 



dlog/X 47T \ 47T 

can be re-written in the form 



dlogfx a s (fi) 

therefore 



47T 47T nn n u! 

2 A) log-, 



and finally a s (/jf) is expressed via a s (/i) as 



<*.M = ^ r ■ (103) 



l + 2A^log^ 

47T /i 



This solution can be written in the form 



2-7T 

a «(/ i ) = 7Z-» ( 104 ) 

A) log 



where Aj^g plays the role of an integration constant (it has dimensionality of energy). 
If higher terms of expansion of (3(a s ) are taken into account, the renormalization group 
equation cannot be solved in elementary functions. 

A surprising thing has happened. The classical QCD Lagrangian with massless quarks 
is characterized by a single dimensionless parameter g, and is scale invariant. In quantum 
theory, QCD has a characteristic energy scale A^jg, and there is no scale invariance — at 
small distances 1/A^jg the interaction is weak, and perturbation theory is applicable; 
the interaction becomes strong at the distances ~ 1/A^jg. This is a consequence of the 
scale anomaly. Hadron masses^ are equal to some dimensionless numbers multiplied by 
Aj^g; calculation of these numbers is a non-perturbative problem, and can only be done 
numerically, on the lattice. 



except the (n^ — l)-plet of pseudoscalar mesons, which are the Goldstone bosons of the spontaneously 
broken SU(nf)A symmetry, and their masses are if n/ quark flavors are massless. 
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Figure 7: a s (/i) from various experiments [13] . 

Values of ct s (/j,) are extracted from many kinds of experiments at various characteristic 
energies /i, see [H]. Their \i dependence agrees with theoretical QCD predictions well 
(Fig. [7]). Of course, all known terms of f3(a s ) (up to 4 loops) are taken into account here^|. 
If these results are reduced to a single /i = mz, they are consistent (Fig. [8]); this fact 
confirms correctness of QCD. 

6 Quark masses 

Until now we considered QCD with massless quarks. With the account of mass, the quark 
field Lagrangian is 

L q = qo(iYD ll -m )q . (105) 

8 Dccoupling effects which arise at transitions from QCD with nf + 1 flavors one of which is heavy to 
the low energy effective theory — QCD with nf light flavors are also taken into account. 
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Figure 8: a s (mz) from various experiments [H]. 

The MS renormalized mass m(/i) is related to the bare one m (appearing in the La- 
grangian) as 

m = Z m {a{n))m{n) , (106) 



(107) 



where Z m is a minimal ( 1731) renormalization constant. 
The quark self energy has two Dirac structures 

S(p)=^(p 2 ) + m S 5 (p 2 ), 

because there is no helicity conservation. The quark propagator has the form 

S(p) 



1 



1 



1 



/-m - /£V(p 2 ) " rn Z s (p 2 ) 1 ~ S y (p 2 ) 1 + S 5 (p 2 



Sy(p 2 



-m 



It should be equal Z q S r (p; //) where the renormalized propagator 5V (p; /x) is finite at e — > 0. 
Therefore the renormalization constants are found from the conditions 



finite 



1 + E 



s 



and hence 



1 + S 5 )Z g Z m = finite . 



finite . 
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At — p 2 ^> m 2 the mass may be neglected while calculating E^Qo 2 ). In the one-loop 
approximation, retaining only the 1/e ultraviolet divergence, we get 

S 5 = C F (3 + a(/i))^. (108) 

Therefore 

OL 

Z m a s ) = 1-3C F — + ■■■ (109) 

47T£ 

The result does not depend on the gauge parameter, this is an important check. 

The dependence of m(/i) is determined by the renormalization group equation. The 
bare mass m does not depend on \i\ differentiating HI 06[) in dlog/x, we obtain 

=^+ 7m (a s (/x))m(/i) = 0, (110) 
a log /i 

where the anomalous dimension is defined as 

d\ogZ m (oL 8 (n)) 
a log// 

For a minimal renormalization constant ( 1731) we obtain from f illip with one-loop accu- 
racy 

lm{a S ) = 7mOT~ H = 2Zl—— H 

47T 47T 

Hence the renormalization constant Z m has the form 

Z m {a s ) = l-— — + ••• 

From (I1U9I) we conclude that 

lm (a s )=6C F -^ + --- (112) 

47T 

Dividing ffTTOl) by flTOOl) (at e = 0) we obtain 

cHogm _ j m (a s ) 
c/loga s 2/3(a s ) ' 

It is easy to express m(//) via m(/x): 

/" Qs ^ 7 m (a s ) da s 
m(//) = m(//) exp / ™ s — - . (113) 

</a s ( M ) 2/3 (a fl ) a s 

Retaining only one- loop terms in 7 m (a s ) and /3(a s ) we get 

m(jjf) = m{p) [^^) ■ (H4) 
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Figure 9: m& from various experiments [15 
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Figure 10: m c (3 GeV) from various experiments [15] . 



26 



Quark masses are extracted from numerous experiments, see the review [15] for b and c. 
For b quark, the quantity fn h is usually presented; it is defined as the root of the equation 

m b (fh b ) = fh b 

(rrib(fj) is the MS mass of b quark), see Fig. [9j For c quark this renormalization scale is 
too low, therefore results for m c (3 GeV) are presented (Fig. [TUl) . 

Acknowledgments. I am grateful to D. Naumov for inviting me to give the lectures 
at the Baikal summer school. 
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